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Abstract 


This is a study to assist in the understanding of earth near 
surface structure. Higher order moments are used to detect the 
density distribution as well as to seek patterns found in 
geological structures. It is shown how higher order moments at 
points outside a mass structure are determined as well as how to 
recover the mass distribution from the higher order moments. It 
is interesting to note that the first moment at a point P outside 
the mass structure,Vo(P), is the entire mass and the second 
moment, Vi(P), is the potential at P due to the mass structure. 
Usually only the mass and the potential function are used to 
determine the density distribution in a body. In this study an 

oo 

infinite function sequence i v n (^n=0* s re Q u i re d to uniquely 
determine the density distribution. 
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Density Determination 

Suppose that A is a material body, i.e., not a point mass, that is 
modeled by assuming that there is a non negative continuous 
function g defined on a closed and bounded subset K of euclidean 
three space E, where K is a geometric approximation of A and g is 
an approximation to the density function of A. Note that K need 
not be connected, i.e., K need not be "in one piece". 

We shall denote the real number line by R, and the euclidean inner 
product on by < • ,• > which induces a norm II • II. If r is a positive 
number and P is a point in E, then 

S r (P) = { Q.is in E: IIP-Q.II<r} 
and BdrySr(P) = { Qis in E: II P - Q. II = r }. 

Since K is bounded there is a number r>0 and a point C, the 
geometric center, of K such that K is contained in S r (C). Let D 
denote the complement of S r (C) in E. 



For each non negative integer n, V n is the real valued function 
defined on D by 


V n (P)-J* E(O) dQ. 
S r (C) II P - Q.H n 


for P in D. 
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The above integral is a triple integral. Notice that Vo(P) is the 
mass of K, and Vi(P) is the Newtonian potential at P due to the 
mass K. 

Suppose g is a non negative continuous function defined on E, Si 
is a ball containing K centered at the geometric center of K, S 2 is a 
ball properly containing Si that is concentric with Si and also g is 
zero on the complement of S 2 . 



There is a one-to-one correspondence between the set 
v={{ V n (P) } n=0 , P is in the complement of sphere S2 } 
and g, moreover g can be constructed from the set V. 

Suppose P is in D. Let mp be the mass function defined on the 
real line R by 

mp(x)= 0 if x < 0 

and mp(x) = / g(Q)dCLifx>0. 

Sx(P) 
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Observe that for each point P in D, mp is a real-valued, non 
decreasing continuous function on R, which has a continuous 
derivative and hence is of bounded variation on R. ( Recall that 
mp is constant except on a finite interval of the line.) 

m'p(x) = 0 if x <0, 

m'p(x) = J* g(x,z) dz if S r (C) intersects BdryS x (P), and 

S r (C) O BdryS x (P) 

m'p(x) = 0 if x> 0 and S r (C) does not intersect BdryS x (P). 

If P is a point, then 

up = inf { II P - QJI : Q.is in S r (C) } - II P - C II - r 
vp = sup { II P - Q. II : Q.is in S r (C) } = II P - C II + r. 

V P 

Observe that V n (P) = J* (l/t n ) dmp(t) for n = 0,1,2, .. . 

Up 
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Making a change of variables s = 1/t for up < t < vp f we have 

Vp" 1 

V n (P) = J* s n dmp(l/s) for n = 0,1,2, . . . 

Up -1 

Define for each point P in D 

Ap= vp-1 , Bp = up'l 

and Mp(x) = M - mp(l/x) Ap<x<Bp. 

Note that Mp(Ap) = 0, and Mp(Bp) = M, 

hence extend Mp to (0, Bp) by defining Mp(x) = 0 for 0 < x < Ap. 
Finally 

M'p(Ap) = m'p(l/Ap)(l/Ap 2 ) = 0 and hence Mp has a continuous 
derivative on (0, Bp). 

We now have 

Bp 

Vn(P) = J* x n dMp(x) for n = 0,1, 2, 3 . . . , 

0 

M’p(Bp) = m' p (1 /B P )( 1 / B P 2) = 0, 

and the sequence {V n (P)ln=0 satisfies the Hausdorf conditions 
(See appendix 1). 

Suppose for each point P on the surface of S 2 , where S 2 is a 
sphere properly containing and concentric with the sphere Si 
which contains K, with center C, (the geometric center of K), the 
sequence 

{ V n (P) } n is known. 

For each such point P consider the line Lp containing P and C the 
center of Si. 
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From the sequence {V n (P)} n we can recover (see Al) the function 
mp. Hence for each point P on S 2 we have mp, then using the 

method of bilinear forms we may recover (see appendix 2) the 
function g. 
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APPENDIX A1 

Suppose that the number sequence [C n ] n= Q is known and satisfies 

the Hausdorf Condition i.e., there is a positive number H such 
that 

n n-t 

X B(n,t) | X B(n-t,s) (-1 )S C t+S I < H for n = 0,1,2, .. . 
t=0 s=0 

where B(n,t) is the binomial coefficient n! /( t! ( n-t )! ). 

For each non negative in n and each number x in the number 
interval [0,1] define 
w n-t 

gn(x) = X B(n,t)X B(n-t,s) (-1 ) s Ct+s if w/n <x< ((w +l)/n). 
t=0 s=0 



I 


The function sequence gn then converges pointwise to a function 
g on [0,1]. The function g is of bounded variation on [0,1] with 
total variation not exceeding H and has the property that 
1 

J* x n dg(x) = C n for n = 0,1,2, ... . 

0 
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APPENDIX A2 

If for each point P on the surface of a hemisphere we know 
Mp(x), then we can determine g on K in the following manner. 

A sequence of functions is created that will converge pointwise to 
g, on a dense subset of S 2 (g is the restriction of g to the dense 
subset ). Since g is continuous, we can extend g to g on S 2 . 

We define for each positive integer n a function fn on the three 
dimensional normalized grid Gn. 
n n n 

f n ( x, y, z) = 2 2 2 G>j(x) Oj(y) Ok(z) wyk for the point (x,y,z) 
i=1 j=1 k=1 

in the grid Gn where n + 1 is the number of grid points on a side 
and h is 1/n. The problem is to determine Wijk for 1 < i,j,k < n. 
Below is a two dimensional schematic to suggest the three 
dimensional case. 



<£ function on grid Typical 0 function 
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The function f n is integrated over the same portion of the grid as 
is intersected by the sphere centered at a point P of radius r (see 
schematic figure below) and set equal to mp(r), which was 
constructed from the sequence {V n (P)}n- By selecting various 
points P and various radii we have a system of linear equations 
in Wjjk that we can solve and thus determine the function f n . As 
stated before, the sequence {f n } n converges pointwise on dense 
subset to g, which is then extended to the function g. 
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